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Abstract 

We study the site-diluted Ising model in two dimensions with 
Monte Carlo simulations. Using finite-size scaling techniques we com- 
pute the critical exponents observing deviations from the pure Ising 
ones. The differences can be explained as the effects of logarithmic 
corrections, without requiring to change the Universality Class. 
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Recently it has been reported using Monte Carlo simulations, that 
the site-diluted Ising model in two dimensions seems to present a second or- 
der transition line with concentration dependent critical indices. Previously, 
other authors [0, had claimed, using analytical methods, an Ising critical 
behavior corrected with logarithms. 

The field theoretical predictions for this model (based on Renormaliza- 
tion Group and Conformal Field Theory) assesses that disorder does not 
change the t] or u exponents since it only changes subleading terms. However, 
the specific heat, at the critical point, diverges as log(logL) in the disordered 
system, while it does as logL in the pure case L being the lattice size. 

In this letter, we extend the methods developed in ref. [Q for the four 
dimensional site-diluted Ising model to the two-dimensional case. 

We observe that, although an apparent variation of the indices seems to 
happen when varying the concentration, this can be explained as a transient 
effect. In fact, a pure Ising value for the indices plus logarithmic corrections 
fits very well our Monte Carlo (MC) data. When preparing this letter, other 
authors have also reported similar conclusions in the bond diluted Ising model 
using different techniques 0. 

The numerical methods as well as the analytical computations based on 
the perturbative renormalization group are very similar in four and two di- 
mensions, so we will report here just the minimal details required, centering 
ourselves in the description of the results. We address to reference 0] for 
further details of the method. 

Our procedure is based on a Finite-Size Scaling (FSS) analysis. We per- 
form MC simulations in the critical region for several values of the concentra- 
tion p. For each p value, we generate hole configurations in a square lattice 
filling the sites randomly with probability p. 

For each hole configuration (sample), we perform a MC simulation of 
the Ising model defined as the set of spins lying in the filled sites coupled 
through a nearest-neighbor interaction. In the smallest concentration we use 
the cluster Swendsen-Wang algorithm ||^ to update the signs of the spins. 
For the rest of concentrations, we found the Wolff single cluster version of 
this algorithm to be more efficient. At each sample, we measure the 
energy, the magnetization and the Fourier transform of the Ising field at 
minimal momenta in about 100 independent spin configurations. We store 
the independent measures in order to compute derivatives with respect to 
the coupling and to extrapolate the results to close values of the coupling 
and the dilution. We have extrapolated in jS in all cases but p ^ 2/3 where 
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L 


p = 1.0 


p = 0.88889 


p = 0.75 


p = 0.6666 


24 


1.009(4) 








32 


1.009(5) 


1.078(6) 


1.142(9) 


1.182(12) 


48 


0.995(4) 


1.072(6) 


1.114(9) 


1.170(12) 


64 


1.001(4) 


1.079(5) 


1.112(9) 


1.151(10) 


96 


1.010(5) 


1.075(6) 


1.111(9) 


1.125(12) 


128 


0.993(4) 


1.066(5) 


1.098(8) 


1.141(13) 


192 




1.065(6) 


1.096(9) 


1.140(13) 


256 


1.004(5) 









Table 1: The v exponent for (L, 2L) pairs at different dilutions. 



a p extrapolation performs better. 

To reduce statistical errors, we have generated 10, 000 hole samples for 
each {L^p) pair, L being the lattice size. We have simulated at concentra- 
tions: p = I, 8/9, 3/5, 2/3. Let us remember that the percolation threshold is 
at Pc — 0.59 0. The FSS method that we use is based upon the ratio of 
several observables: magnetization, susceptibility, correlation length. Binder 
cumulants, and their derivatives for two different lattice sizes Li and L2. At 
the parameter values where the correlation lengths ratio matches L2/L1 we 
expect that, in the absence of logarithmic corrections 

nn'^i'^l = iL,/L,ro/u + 0(Lr ) , (1) 
0(Li,/3,p) 

where xq is the critical exponent for the observable O, e.g. 7 for the sus- 
ceptibility, u for the correlation length, etc. uj is the universal corrections- 
to-scaling exponent. The /9-derivatives of the previous observables go with 
the corresponding exponents plus l/u. In all cases we use pairs of lattices of 
sizes L and 2L. 

In table 1 and 2 we report the results for the exponents u and t] respec- 
tively, using the relation (|l]) for the /?-derivative of the correlation length in 
the former case and the susceptibility in the latter (after applying the scaling 
relation 77 = 2 — 7/1/). Notice that there is some statistical anticorrelation 
between sucessive even (odd) rows. We recall the critical exponents in the 
/3 — > 00 limit (pure site percolation), u = 4/3 and t] = 5/24 conjectured by 
Nienhuis and in agreement with a recent MC work [jl2| . 

Although the u exponent seems to be non-constant as a function of the 
concentration, the values for rj are very stable. This fact was also observed 
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p = 0.75 


p = 0.6666 
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0.2465(7) 








32 


0.2466(8) 


0.2495(9) 


0.2504(10) 


0.2454(18) 


48 


0.2495(8) 


0.2487(7) 


0.2498(11) 


0.2448(14) 


64 


0.2499(8) 


0.2490(8) 


0.2469(8) 


0.2455(15) 


96 


0.2498(8) 


0.2501(8) 


0.2460(9) 


0.2456(14) 


128 


0.2497(8) 


0.2495(8) 


0.2469(9) 


0.2461(13) 


192 




0.2485(8) 


0.2483(7) 


0.2483(14) 


256 


0.2517(7) 









Table 2: The t] exponent for (L, 2L) pairs at different dilutions. 



Pc 


/3c 


UJ 


xVd.o.f. 


Li - {L2 - L2) 


1 


0.440682(5) 


1.5(8) 


15.2/12 


24-(32-512) 


0.88889 


0.53781(2) 


1.0(3) 


2.6/8 


48-(64-384) 


0.75 


0.77125(8) 


0.9(5) 


7.7/6 


48-(96-384) 


0.66661(3) 


1.10 


0.6(2) 


9.6/8 


32-(64-384) 



Table 3: Critical parameters for several dilutions. We also show the computed 
corrections to scaling exponent u, the fit quality and the range of lattice sizes 
used (crossings from Li — L\ to Li — Lg). 



in the scaling of the Yang-Lee zeroes at the critical point |T3| 



An interpretation of these results could be a continuous set of Universality 
Classes, as assumed in reference [Q. However, a simpler scenario is a fixed 
value for v (that of the Ising Model) plus logarithmic corrections. 

Using the results from refs. p|, |14[ it is possible to show that the derivative 
of the correlation length with /3 reads at critical point as 

^^^"-Vl+l-logL0 + 2 + 2ClogL ) 

where C depends on the dilution. Computing the v exponent from (|1]) and 
(H) we obtain an apparent behavior as 



V 



.apparent _ (3) 



In figure 1 we plot the values of v reported in table 2 as a function of 
1/ logL. All our MC data are compatible with a behavior of type (|]). 
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Figure 1: v exponent as a function of 1/ logL for several values of p. 



The critical parameters can be computed with great accuracy studying 
the crossing of several quantities (second momentum correlation length di- 
vided by the lattice length, and Binder parameter for the magnetization) for 
the different lattice sizes. The deviation of the crossing point for lattices sizes 
L, sL scales as 

1 — 

Ap^ oc ^^^^L— (4) 



We address to reference 10 for details of the method. The results are re- 



ported in table 1^. Notice that for the Ising Model the method gives the correct 
value within errors for the critical temperature (1 part in 10^). Regarding 
the oj exponent, whose conjectured value is u; = 4/3 [|1^] for the pure system, 
to obtain more accurate results, simulations on smaller lattices should be 
added, but it would increase the systematic errors in the determination of 
the critical coupling. 

We finally measure the specific heat Cy at the critical points. The results 
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Figure 2: Specific heat as a function of the double logarithm of the lattice 
size. The lines correspond to linear fits for L > 96. The three lines for 
each case show the effect of the error in the determination of the critical 
parameters. 

are displayed in figure ^ Excluding the Ising limit, a linear function of 
log(logL), as predicted in ref. 0, fits very well the data. 

We thus conclude that the site diluted Ising model in two dimensions 
belongs to the same Universality Class as the pure Ising model, although we 
find strong logarithmic effects. 

We thank to the CICyT (contracts AEN94-0218, AEN96-1634) for partial 
financial support. JJRL is granted by EC HMC (ERBFMBICT950429). 
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